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ABSTRACT 


Digital  simulation  and  modeling  of  a  turbine  engine  horsepower 
extraction  system  is  presented.  The  system  consists  of  a  digital 
valve,  a  fluid  transmission  line  (piping),  and  a  waterbrake.  The 
system  model  was  developed  using  the  basic  lumped-parameter  approach. 
A  computer  modeling  language  called  Advanced  Continuous  Simulation 
Language  (ACSL)  represents  the  model's  response  to  step  and  profile 
Inputs.  The  simulation  of  step  and  profile  inputs  demonstrate  the  con¬ 
trol  system's  ability  to  set  power  levels  and  control  changes  much 
faster  and  with  better  accuracy  than  a  manual  system.  The  waterbrake 
equation  was  empirically  derived  and  validated. 
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Introduction 


Modern  control  systems  are  necessary  to  control  complex  tasks 
effectively.  A  human  operator  is  unable  to  correct  for  error  in  milli¬ 
seconds.  A  digital  computer  can  accomplish  complex  tasks  in  milli¬ 
seconds.  Many  of  today's  computer  controllers  provide  reliable  and 
consistent  error  corrections  on  both  complex  and  repetitive  tasks.  To 
control  a  system  effectively,  whether  a  heating  system  or  rocket  guid¬ 
ance  system,  the  system  must  be  understood.  Modeling  subdivides  a 
complex  system  into  its  basic  components.  Once  a  model  is  developed, 
whether  digital  or  analog,  one  needs  to  evaluate  the  control's  ability 
to  regulate  the  system. 

This  thesis  models  a  gas  turbine  horsepower  extraction  system 
(waterbrake)  using  a  computer  modeling  language  called  Advanced  Con¬ 
tinuous  Simulation  Language  (ACSL).  The  control  requirements  for  a 
waterbrake  system  are  stated  in  Chapter  I.  Chapter  II  develops  the 
horsepower  extraction  (HPX)  model  from  the  basic  parameters  and  system 
equations  for  each  component.  The  computer  program,  written  in  ACSL, 
solves  the  system  equations  in  Chapter  III.  The  model  response  to 
various  inputs  -  step  and  profile  inputs  are  shown  in  Chapter  IV. 
Chapter  V  summarizes  the  results  and  provides  a  recommendation. 
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CHAPTER  I 


CONTROL  REQUIREMENTS 

Background 

Most  of  the  turbine  engine  horsepower  extraction  (HPX)  systems  at 
AEDC  use  a  waterbrake  connected  to  the  engine  gearbox  shaft  to 
simulate  aircraft  accessory  loads.  Presently  the  HPX  Is  set  manually 
which  Is  time  consuming  and  does  not  provide  accurate  simulation  of 
aircraft  accessory  loads  during  transient  engine  operations.  Setting 
the  HPX  automatically  would  save  a  few  minutes  during  each  airstart 
test.  Testing  turbine  engines  at  Arnold  Engineering  Development  Center 
(AEDC)  Is  expensive  and  a  few  minutes  saved  each  test  applied  to  the 
total  number  of  tests  means  a  potential  saving  (for  thirty  tests)  of 
approximately  $35, 000. 00/year.  Presently,  the  HPX  1s  set  at  one  engine 
shaft  speed  (rpm).  This  setting  Is  not  adjusted  for  changes  In  shaft 
speed.  As  a  result,  the  horsepower  extracted  does  not  match  the  typi¬ 
cal  aircraft  accessory  loads.  A  computer  controller  can  automatically 
correct  for  changes  in  shaft  speed  and  adjust  water  flow  to  the  water- 
brake  during  engine  simulated  airstarts.  Setting  the  HPX  level  cor¬ 
rectly  and  continuously  would  significantly  reduce  the  difference 
between  test  simulated  HPX  and  actual  aircraft  HPX  (Figure  1). 

Waterbrake 

The  waterbrake  Is  a  key  component  in  the  HPX  system.  The 
waterbrake  1s  a  power  absorption  device  driven  by  a  shaft  from  the 
engine's  gearbox  (Figure  2).  Power  1s  absorbed  by  water  as  It  passes 
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2  4  6  8  10  12 

Shaft  Speed  x  1000  (RPM) 

Figure  1.  Typical  Aircraft  Horsepower  Loads  Versus  Simulated 
Loads 
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Figure  2.  Waterbrake  System  Hardware 


through  a  series  of  rotating  discs  and  stationary  plates.  The  water  is 
discharged  from  the  waterbrake  at  a  higher  temperature  due  to  the  basic 
energy  absorption  process.  The  level  of  power  absorbed  by  the  water- 
brake  can  be  regulated  by  controlling  the  water  flow  rate  or  shaft 
speed  [1]^. 

The  waterbrake  horsepower  extraction  system  uses  reaction  torque 
and  rotor  speed  measurements  to  determine  horsepower.  Measuring  the 
torque  and  the  speed  provides  almost  Instantaneous  readout  capability. 
The  waterbrake  system  hardware  shown  in  Figure  2  uses  a  trunion 
bearing,  a  lever  arm,  and  a  load  cell  to  measure  torque.  The  engine's 
gearbox  shaft  speed  is  determined  from  the  engine  core  speed  pickup  and 
a  known  gearbox  ratio.  The  correct  size  of  waterbrake  required  is 
determined  from  the  brake's  power  absorption  capacity  curves.  Figure  3 
shows  the  maximum  HPX  versus  a  range  of  shaft  speeds  (rpm)  for  typical 
waterbrakes  used  at  AEDC. 

System  Operation 

The  current  horsepower  extraction  system  is  manually  operated  as 
shown  on  the  top  of  Figure  4.  A  one-inch  water  supply  line  feeds 
water  at  65  psig,  into  a  control  valve.  A  0.5-inch  pipe  transports  the 
water  to  the  waterbrake  15.5  feet  away;  where  it  passes  through  the 
waterbrake  and  is  drained.  A  load  cell  measures  the  reaction  force 
using  a  trunion  bearing  and  linkages  (Figure  2).  The  extracted 

lumbers  in  brackets  refer  to  similarly  numbered  references  in  the 
List  of  References. 
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Figure  3.  Waterbrake  (Hydra-brake)  Horsepower  Curves 
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a.  Manual  controlled  system. 


Figure  4.  Comparison  Between  Manual  and  Computer  Controlled  Systems 


horsepower  is  calculated  from  the  reaction  force,  linkage  length,  and 
gearbox  shaft  speed.  The  automatic  control  system  on  the  bottom  of 
Figure  4  constantly  corrects  to  the  desired  horsepower  by  using  a  com¬ 
puter  program  instead  of  an  operator  to  calculate  power  and  automati¬ 
cally  adjust  water  flow. 
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CHAPTER  II 


MODELING  THE  HORSEPOWER  EXTRACTION  SYSTEM 


Control  System 

The  primary  control  system  goal  is  to  set  the  desired  horsepower 
accurately  during  steady-state  engine  operations  which  saves  expensive 
engine  testing  time.  The  secondary  goal  is  to  provide  relative  control 
during  an  engine  air  start  transient  characterized  by  rapid  rotor  speed 
increase  and  step  increases  in  horsepower  demands.  A  mathematical 
model  was  written  representing  the  physical  system  analytically,  using 
a  series  of  equations.  A  computer  model  simulation  will  help  determine 
if  the  physical  system  can  meet  the  primary  and  secondary  goals 
mentioned  above. 

The  horsepower  extraction  (HPX)  control  system  has  four  -  main 
components:  control  valve,  transmission  line,  waterbrake  and  con¬ 
troller.  Each  one  of  the  main  components  can  in  turn  be  represented  by 
a  lumped-parameter  model. 

Modeling  Technique 

The  lumped-parameter  technique  assumes  the  physical  distributed 
characteristics  of  pipe  wall  friction  (resistance,  RE),  fluid  com¬ 
pressibility  (capacitance,  CE),  and  water  momentum  (inductance,  LE). 
along  the  pipe  can  be  estimated  by  one  lump  (value)  for  each  of  the 
three  characteristics  (RE,  CE,  LE).  By  reducing  the  distributed 
characteristics  of  a  pipe  into  lumped-parameters  (RE,  CE,  LE)  the 
complex  and  cumbersome  equations  of  the  flow  can  be  simplified. 
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Further  information  on  this  topic  can  be  found  in  Reference  12].  The 
lumped-parameter  approach  reduces  the  fluid  (physical)  system  into  its 
basic  fluid  parameters.  The 
correlation  between  the  fluid  terms  and  electrical  terms  are  shown  in 
Table  1.  The  three  basic  fluid  parameters,  equivalent  resistance 
(RE),  equivalent  capacitance  (CE),  and  equivalent  inductance  (LE),  are 
determined  for  the  valve,  transmission  line,  and  the  waterbrake. 

The  fluid  equations  and  basic  parameters  are  developed  for  the 
transmission  line  on  the  next  few  pages. 


Equivalent  Resistance  fREI 


Calculation  of  equivalent  resistance  (RE)  shown  in  Figure  5  is 
computed  from  Bernoulli's  Equation  [3], 


where 

H  =  head  loss  (ft  Ibf/lbm)  0  =  pipe  diameter  in  inches 

f  =  friction  factor  v  =  fluid  velocity  velocity  in  ft/sec 

L  =  pipe  length  in  inches  gc  =  32.2  ft  Ibm/lbf  sec2 

p  =  density  in  Ibm/ft^  dp  =  change  in  pressure  (Ibf/ft^) 
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Table  1.  Correlation  of  Mechanical  and  Electrical  Terms 


Terms 

Electrical 

Fluid  Equivalent 

Symbols 

Description 

Resistance 

R 

RE 

fL/(D  2g  p  A2) 

Capacitance 

C 

CE 

pV/B 

Inductance 

L 

LE 

L/A  gc 

Voltage/Pressure 

V 

P 

pressure  (p) 

Current/Mass  Flow 

i 

w 

Mass  Flow  (w) 

Voltage  Drop/ 
Pressure  Drop 

dv 

dp 

pressure  drop 

where 


L  =  pipe  length 
D  =  pipe  diameter 
p  =  water  density 
A  =  pipe's  cross-sectional  area 
gc  =  conversion  factor 
f  =  frictional  factor 
B  =  bulk  modulus 
V  =  vol ume 
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t  Fluid  Transmission  Line 

P,  LE  RE  P2 


III  i 


Equivalent  Electrical  Transmission  Line 

PiandPj  =  Ibf/in2  Analogous  to  Voltage 
WjandWj  =  Ib^/sec  Analogous  to  Current 
Wi  =  W|,  =  Wj  =  Flow  Through  Inductor  and  Resistor 
Wj  =  Wj  -  W;  =  Flow  Through  Capacitor 


Figure  5.  Fluid  to  Electrical  Transmission  Line  Analogy 


(2) 


H  = 


P 


and 


However, 


and 


-P,  =  Pf 


f  W  \ 

u  =  —  j  ,  wherew  =  flow  mtelbm/sec. 


P^Ahgc 


(3) 

(4) 

(5) 


If  we  let 


RE  = 


C-) 


pAhgc 


(6) 


The  equivalent  resistance  (RE),  of  the  transmission  line  Is 
calculated  below. 


(7) 


where 

L2  =  effective  length  of  transmission  line  pipe  =  15.5  feet 
D2  =  Internal  pipe  diameter  =  0.5  Inch 
Pw  =  water  density  =  62.4  Ibra/ft^ 

A2  =  cross-sectional  area  of  the  pipe  =  022/4 
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A2  =  0.196  Sq  inches  =1.36  E-3  Sq  feet 
Qj.  =  gravitational  constant  =  32.2  Ibm  ft/lbf  sec 
f  =  frictional  factor. 

To  determine  the  frictional  factor  (f)  one  needs  the  Reynolds 
Number  (RN)  and  the  pipe  roughness  e/D  [4]. 

e  =  roughness,  commercial  =  0.00015 

e/D2  =  0.00003 

RN  =  Pw  V2  D2/(pw  gc)  (8) 

pw  =  absolute  viscosity  of  water  @  T  =  68*F 
=  2.09  E-5  Ibf  sec  /sq  ft 

V2  =  average  velocity  =  average  volume  flow/A2  (9) 

=  8  gal/min  (ft3/7.48  gal)  (min/60  sec)/1.36  E-3  ft 
=  12.96  ft/sec 

RN  =  4.98  E  4  (turbulent  flow) 
f  =  0.022  from  Moody  chart  [3]. 

Knowing  the  frictional  factor  we  can  now  calculate  the  equivalent 
resistance  (RE2), 

RE2  =  7.61  (Ibf/in2)/(lbm2/sec2). 

Equivalent  Capacitance  (CEl 

The  next  constant  to  be  determined  is  the  capacitance.  The  line 
capacitance  is  attributed  to  the  compression  of  the  fluid  (water),  com¬ 
pression  of  the  entrained  air  and  the  expansion  of  the  pipe.  Since  the 
pressure  in  the  pipe  is  very  low  (less  than  65  psi),  the  pipe  expansion 
can  be  neglected.  The  amount  of  air  entrained  in  the  water  is  not 
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known  but  considered  to  be  between  0  to  a  few  percent.  Some  air 
entrainment  Is  typical  In  a  water  system  and  one  percent  by  volume  was 
used  In  the  simulation  model.  The  equivalent  capacitance  (CE)  In 
Figure  5,  page  12,  can  be  computed  by  considering  the  Bulk  Modulus  (B). 


CE^  = 


Pu,^2!^2 

B 


(10) 


where  B  =  bulk  modulus  water  with  entrained  air  one  percent  by  volume 
1/B  =  1/Bw  +L/Ba  [2] 

1/B  =  (0.99/3.25E5  +  0.01/65) 

B  =  6374  Ibf/sq  in. 

Therefore, 

CE2  =  2.13  E-4  Ibm  sq  in/lbf 


The  equivalent  capacitance  (CE)  Is  derived -below: 


and 


therefore. 


dp 

B  =  —V 
dv 


V  =  volume  ft 


dV=&V 

B 


B 

dp  =  -dV 


^  _  B  ^ 

dt~  V  dt 


(11) 


(12) 
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However 


dV  li) 
dt  p 


Defining 


then 


—  =  —  li; 

dt  Vp 


=  h  \o 


(13) 

(14) 

(15) 


Equivalent  Inductance  (LE21 

Inductance  (Figure  5,  page  12)  can  be  calculated  by  considering 
Newton's  second  law  [2,  5],  and  assuming  RE  =  0,  then 


F  = 


mo 


where, 

gc  =  32.2  ft  Ibm/lbf  sec2. 
F  *  force,  Ibf, 


m  =  mass  of  fluid  Ibm, 
a  =  acceleration  of  the  fluid  in  ft/sec2 


then 


F  =  APA 


where 

£J>  =  P,-P,. 

The  mass  of  fluid  (m)  in  the  line  of  length  L  is. 


(16) 


(17) 
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m  =  pAL. 


(18) 


Let 


then 


Define 


then 


therefore 


Let 


,  dv 

accelerationa  =  — 
dt 


ma  pAL  dv 

F=  —  =  APA  =  - - 

8,  8^  dt 


_  w  ^  dv _ ^  d^ 

pA  dt  pA  dt 


APA  = 


pAL  1  dw 
pA  dt 


L  dw 

AP= - . 

g^A  dt 


where  L2  =  pipe's  length  =  15.5  feet 
gc  =  conversion  factor 
A2  =  pipe's  cross-sectional  area 
=  0.196  square  inches. 


(19) 


(20) 


(21) 


(22) 


(23) 
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Fluid  Transmission  Line 


The  fluid  transmission  line  consists  of  five  sections  of  pipe,  and 
three  90  elbows  which  transport  the  water  from  the  digital  valve  to  the 
waterbrake.  The  fluid  to  electrical  transmission  line  analogies  are 
presented  in  Figure  5,  page  12,  and  the  description  of  the  equations 
are  shown  below. 

The  fluid  transmission  line  is  converted  to  a  electrical  circuit 
to  model  a  control  system.  The  same  techniques  used  to  solve  an  elect¬ 
rical  circuit  will  apply  to  the  analogous  fluid  circuit.  The  pressure 
drop  (analogous  to  voltage)  drop  from  Pi  to  Pz  Is  equal  to  pressure 
drop  across  the  Inductor  (LE)  plus  the  pressure  drop  across  the 
resistor  (RE)  (2,5,6]. 

P^-P^  =  LE^  ~  (24) 

Electrically  this  would  be  analogous  to 

=  (25) 

solving  for  W2 

The  pressure  P2  is  defined  by  the  Equation  (28)  using  the  fact 
from  the  law  of  conservation  of  matter. 
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Water  stored  equals  water  in  minus  water  out  is  shown  in  Equation  (27) 
below. 


dP 


Rearranging  terms  to  solve  for  P2 


(27) 


(28) 


Waterbrake  Power  Equation 

The  waterbrake  is  a  compact  power  absorption  device  used  to 
simulate  loads  on  turbines  and  other  rotating  equipment.  Equation  (32) 
below  calculates  the  horsepower  absorbed  assuming  the  power  is  dissi¬ 
pated  by  acceleration  and  deceleration  of  the  fluid.  Figure  6  shows  a 
cross-sectional  view  of  the  waterbrake  with  the  radius  to  the  shaft, 
gap,  and  holes. 


Power  =  Torque  x  Angular  Velocity  =  T  x  6  (29) 

and 

Torque  ,=  force  x  distance  =  F  x  r  (30) 

since 

Force  =  mass  x  acceleration  =  w  dv/dt  *  w  6  r  (31) 

gives  us 

T  =  CFw  w  er2 

The  total  torque  equals  rotor  torque  from  the  water  passing  (wh) 
through  the  holes  (rh)  plus  water  passing  (wg)  through  the  gap  (rg) 
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Figure  6.  Cross-Sectional  View  of  the  Waterbrake 


minus  initial  torque  from  the  water  passing  over  rotor  shaft  (rs). 
This  equation  is  written  out  below  for  six  stages: 

T  =  CFw  {wh  rh  (rh  6)  +  Wg  rg  (rg  6) 

-  ws  rs  X  (rs  6)}  x  2/stage  (6  stages)  (32) 

where 

rh  =  radius  from  shaft  centerline  to  the  holes  (2.50  inches) 
rg  =  radius  from  shaft  centerline  to  the  gap  (3.00  inches) 
rs  =  radius  of  the  shaft  (1.00  inches) 

stage  =  rotor  +  stator  (shown  in  Figure  6) 

Q  =  rotor  speed  (N)  x  1.144  (gearbox  multiplier),  rpm 
w  =  water  flow  in  gal/min 
wt  =  Ws  =  total  water  flow 
Ws  =  water  flow  around  the  shaft  at  radius  rs 
Wh  =  water  flow  through  the  holes  at  radius  rh 
wg  =  water  flow  through  the  gaps  at  radius  rg 

The  total  flow  through  the  waterbrake  is  equal  to  the  flow  through  the 
hole  plus  the  flow  through  the  gap  as  written  in  Equation  (33) 

Wj  =  Wh  +  Wg  (33) 

Let  Wt  =  >«  w  +  y  w 

where  x  =  percentage  of  total  water  through  the  rotor  holes 
y  =  percentage  of  total  water  through  the  rotor  gap 
Based  on  flow  area  of  the  rotor  holes  and  the  circumferential  gap, 
and  the  equation  fit  to  the  waterbrake  data.  The  x  and  y  values  were 
chosen  to  be:  x  =  0.20  and  y  »  0.80. 
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Factoring  out  6,  gives  the  torque  equation 


T  =  CFw  w  6(xrh2  +  yrg2  -  rs2)  (34) 

The  waterbrake  horsepower  can  now  be  determined  from  the  power 
equation  shown  below. 

Power  =  [CFp  6^  w  (x  rh^  +  y  rg^  -  rs^)]  6  stages  (35) 

CFp  =  conversion  factor  =  5.99  E-10. 

Table  2  compares  the  experimental  test  data  to  the  analytical  HPX 
values  for  various  water  flow  rates  and  rotor  speeds.  The  comparison 
is  incomplete,  due  to  the  lack  of  test  data. 

Waterbrake  Flow  Rate  Equation 

The  flow  rate  (W)  is  very  complex  and  not  easily  modeled.  For 
this  reason,  an  empirical  equation  was  derived  from  the  test  data 
available.  The  flow  equation  is  as  follows: 

w  =  (dP/RE3)l/2  (36) 

where 

w  =  water  flow  through  the  waterbrake, 
dP  =  pressure  across  the  waterbrake, 

RE3  =  resistance  which  is  a  function  of  both  inlet  pressure, 
rotor  speed  and  flow  rate. 

Data  used  to  determine  RE3  can  be  found  in  Figure  7.  The  induc¬ 
tance  (LE3)  and  capacitance  (CE3)  waterbrake  parameters  were  not  con¬ 
sidered  due  to  lack  of  dynamic  test  data  to  characterize  waterbrake. 
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Table  2.  Actual  Waterbrake  Data  Versus  Analytical  Horsepower  Values 


Rotor  Speed 
(M) 
rpm 

Flow  Rate, 

Ibm/sec 

Equation 

Output, 

HPX 

hp 

Experimental 
Data,  HPX 
hp 

Different 
Equation- 
Exp.  HPX 
hp 

2.22 

189 

190 

-1 

12.500 

2.09 

164 

180 

-16 

11,800 

0.76 

54 

80 

-26 

11,800 

1.25 

88 

110 

-22 

I  tjjjllll 

0.70 

39 

45 

-6 

2,500 

1.95 

6 

10 

-4 
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Pressure 
Drop  Across 
Waterbrake 
(psig) 


Waterbrake  Flow  ( Ibm  /sec) 


Figure  7.  Empirical  Waterbrake  Data  with  Lines  of  Constant 
Resistance  at  Designated  Rotor  Speeds 
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Digital  Control  Valve 

A  digital  twelve  element  control  valve,  manufactured  by  The  Digi¬ 
tal  Control  Company  [7],  was  chosen  to  regulate  the  water  flow.  The 
high-tech  valve  has  a  (full  travel)  response  time  of  less  than  100 
milliseconds,  no  hystersis,  and  a  linear.CV  trim.  The  general  form  of 
the  control  valve  equation  is  shown  below  [6]. 

li;  =CF^Cv(Vdp)  (37) 

where 

dP  =  pressure  across  the  valve  (psig) 

CV  =  flow  coefficient  defined  as  the  volume  of  water  In  gallon/ 
minute  that  will  flow  with  a  one  psi  pressure  drop 

CFv  =  conversion  factor. 

The  digital  valve  was  sized  for  a  maximum  water  flow  of  16  gpm,  a 
9  psig  across  the  valve,  and  a  flow  resolution  of  0.015  gpm.  The 
twelve  elements  (solenoids)  each  have  separate  orifices  that  act  as 
switches  when  powered  by  100  psig. 

The  capacitance  for  the  1.0-inch  valve  that  is  12  inches  long  can 
be  found  in  Table  3.  The  valve  inductance  of  (0.04)  is  neglected 
compared  to  the  transmission  line  inductance  of  (2.46),  therefore,  and 
will  be  ignored  to  reduce  the  complexity  of  the  control  system  simula¬ 
tion. 
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Table  3.  Summary  of  Calculated  Parameters 


Section  A;  Valve 

REl  =  f  (valve  area,  dp) 

LEi  =  0.040  (lbf/in.2)/(lbm/sec2) 
CEl  =  5.23  E-5  (lbm)/(lbf/1n.2) 

Section  B:  Transmission  Line 

RE2  =  7.61  (Ibf/1n.2)/(lbm2/sec2) 
LE2  =  2.46  (lbf/in.2)/(lbm/sec2) 
CE2  =  2.13  E-4  (lbm)/(lbf/1n.2) 

Section  C;  Waterbrake 
RE3  =  f(dp,  6,  w) 

Not  considered 
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System  Equations 

The  system  equations  and  a  flow  schematic  of  the  digital  valve, 
transmission  line  and  waterbrake  are  shown  on  Figure  8.  The  electrical 
electrical  analogy  of  the  mechanical  system  Is  also  shown  on  Figure  8, 
Summary  of  the  calculated  parameters  (CE,  RE,  and  LE)  are  provided  in 
Table  3.  Equations  (7),  (10),  and  (23)  define  the  mechanical  equ1v> 
alent  resistance  (RE),  capacitance  (CE),  and  Inductance  (LE).  The 
transmission  line  Equations  (26)  and  (28)  define  the  flow  (w)  and 
pressure  (P)  terms.  The  waterbrake  flow  rate  (w3)  Equation  (36)  was 
derived  empirically  from  the  experimental  data.  The  digital  valve 
water  flow  (wi)  and  downstream  pressure  (P2)  were  determined  by  Equa¬ 
tions  (26)  and  (28),  respectively. 
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/t(W2  -W3)  dt 
0 


W3  =  f  (N.P2.P3.  RKa) 


NOTE:  REj.LEj.CEjandCE,  were  determined  by  equations  7, 23  and  10. 
The  Equations  for  W, ,  and  W,  are  equations  37, 26  and  36. 

The  Equations  P,  and  Pj  are  both  from  equation  28. 


Figure  8.  Horsepower  Extraction  (HPX)  Parameters,  Locations,  and  System  Equations 


CHAPTER  III 


COMPUTER  SIMULATION  PROGRAM 

Background 

The  system  equations  can  be  solved  by  several  approaches.  Some 
approaches  are:  analog  computer  simulation,  custom  numerical  computer 
simulation,  or  Advance  Continuous  Simulation  Language  (ACSL).  The 
extreme  flexibility  and  computing  power  of  ACSL  make  it  the  most 
attractive  approach.  The  highlights  of  the  program  are:  uses 
conventional  Fortran  statements,  has  nonlinear  function  generation  up 
to  three  variables  and  flexibility  to  plot  model  behavior  with  various 
forcing  functions.  The  ACSL  manual  [81  provides  the  detail  of  its 
capabilities. 

Computer  Program 

The  system  equations,  shown  in  Figure  8  were  programed  using  ACSL. 
The  computer  program  is  given  in  Table  4.  The  program  is  separated  into 
four  sections:  valve,  transmission  line,  waterbrake  and  controller. 
The  flow  chart  (Figure  9)  and  system  equations  (Figure  8)  show  the 
interdependence  of  each  component.  The  beginning  of  the  program  pro¬ 
vides  the  initial  conditions,  the  parameter  constants  (CE,  RE,  LE),  and 
the  look-up  tables  of  nonlinear  data  (shown  graphically  in  Figures  10- 
13).  Figure  11  shows  at  lower  rotor  speeds  a  greater  change  in  water 
flow  is  required  to  obtain  the  same  change  in  the  amount  of  horsepower. 
The  physical  system  is  bound  between  0.05-65  psi  and  0.01  (leakage)  to 
2.30  Ibm/sec.  The  valve  position  range  is  0.1  to  98  percent  of  full 
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Table  4.  Computer  Program  Listing 


•••••••••advanced  continuous  simulation  lanquage******** 

ACSL  TMANSLATOR  VEKSZON  4  LEVEL  SL1  SO/1S4  IS. 01. 02  I 

■nST  CASE  WITH  CONVERSION  TO  ASCL" 

■  ■ 

CINTERVAL  CINT  b  .1 

MAXTERVAL  NAXT  •  .010 
■  ■ 

NSTEPS  MSTP  B  1 

PROGRAM  HATERMIAKE  (HPX)  CONTROLLER 
H  a 

-TABLE  OF  VALUE  POSITION  (INPUT)  VS  VALVE  FLUID  FLOW,  PVI(POSIT)  « 
TABLE  PYI.1,6 

/O.O. 10.0.30. .60. .80.0.200. 

.0.0.. 10.. 30. .60.. BOO. 2. 00  / 

a  a 

a  a 

■TABLE  OF  ROTOR  SPEED  (N2,RPMS)  VS  CHANGE  HPX,  FV2(CHPX}- 
TABLE  FY2.1.B 

/3000. ,BBOO. ,8200. , 10000. , 12600. 

.75. .38. .18. .12. .8.600  / 

a  a 

a  a 

■TABLE  OF  ROTOR  SPEED  (N2.RPMS)  VS  HATERBRAKE  RESISTANCE  FV3(R4)- 

n  a 

TABLE  FY3,1,5 

/1000. .2500.0, 10000. , 11800. , 13000. 

.8. 1.6. 3, 4. 0,3.0  ,2.9  / 

■SET  CONSTANTS  ■ 

■INITIAL  PRESSURES  &  FLOWS  ” 

■TABLE  OF  ROTOR  SPEED  (N2,RPMS)  VS  DESIRE  HPX  FV4(HPX0}* 

M  a 

TABLE  FV4,1,9 

/3500. ,4000.0,4010. ,5000. ,5010. ,6000. .6010. ,7000. , 10000.  . . . 

, 10. , 10. ,20. .20. ,30. ,30. ,50. .50. . 100.  / 

■SET  CONSTANT- 

CONSTANT  PSa60.00,P1lB.01.P2I-.01 
CONSTANT  HOI-0.00,W1Ia0.01,H2lB.010 

a  a 

■INITIALIZE  HPX  VALUES  A  VALVE  POSITION  ■ 

CONSTANT  X0L0«.01,CVb.01.TRS>1.0,CHb3.BBE-7  .OC-1.144 
,XDb.01,T0L>  .100 

■  SET  RESISTANCE  (RE),  INDUCTANCE  (LE)  5  CAPACITANCE  (CE)  ■ 
CONSTANT  RE2a7.B10,RBa20.00 

CONSTANT  LE2a2.46 

CONSTANT  CE2a2.13E-4,  CE1aB.23E-B 

a  a 

■  BOUND  THE  SYSTEM  BY  PHYSICAL  LIMITATIONS  " 

CONSTANT  PBb.B0,PTb65.0,P01Bb.20.PD1Tb  60.  ... 

,PDB—B0.0,P0T-80.0  ,C6b-B00.  .CTbBOO.  ... 

,VBa-1.,  VTb  88.  .NBb.00B.HTb2.30 

a  a 

DYNAMIC 

DERIVATIVE 

a  a 

PS-  65. 

■N2a10300.  * 
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Table  4.  Concluded 


•••••••••advanced  continuous  simulation  LANQUAAE*' 

ACSL  TRANSLATOR  VERSION  4  LEVEL  SLI  80/134  12.01.02 

N2>  4000. 4  380.*IUMP(TRS) 

■  ■ 

CV>  FVKPOSIT)  *1.0 
CR«  FV3(N2)  *1.0 
POD-  PS-P1 

■■  SECTION  A  :  VALVE  ■ 

P10-B0UND( P01B , PD1T , POD) 

Ml 1 - . 85«CV*SQRT ( ABS( P 10 ) ) 

W1  ■  B0UN0(WB,WT,W11) 

■  M 

DP12-(W1-W2)/CE1 

OP 1 -BOUNO( PDB , POT , DP 12 ) 

P1-L1NINT(DP1,P1I,PB.PT) 

N  11 

■  SECTION  B  :  TRANHISSION  LINE  " 

C8-  PV2(N2)  *1.0 

CBb  FVKPOSIT)  *1.0 

OWE- ( ( P 1 -P2 ) - ( H2«ABS ( W2 ) •RE3 ) )/LE2 

H2-LIMINT(0W2,W21,WB,WT) 

■  ■ 

0P23-(W2-W3)/CE2 

0P2 -BOUND ( POB , POT , 0P23 ) 

P2-LIM1NT(DP2,P2I.PB,PT) 

■  m 

P2 1  -  ( P2-f  ABS(  P2  )-«■ .  0001 ) /2 . 0 
■SECTION  C  :  NATERBRAKE  ■ 

W33-  SQRT(ABS((P21)/Fy3(N2))) 

W3-  BOUNOCWB.MT.MSS) 

■  ■ 

■SECTION  0  :  CONTROLLER  ■ 

WC>(W3-fABS(W3))/2. 

HPXA- ( N2«aC ) • *3 . •HC»CW 
■HPXD>  10.0  *  20.*RAMP(TRS)  ■ 

HPXE-HPXD-HPXA 

■HPXD*  100.  " 

t»XD>  FV4(N2) 

CHPX-  FY2(N2)  *1.0 
DWV-HPXE«CHPX  «.400 
CP-  B0UND(CB,CT,0HV) 

»  ■ 

■  • 

POS-  LIMINT(CP.XOLD.VB.VT) 

POSIT-  OELAY(P0S,XO,TDL.100O) 

END  $ 

CONSTANT  TSTP  -30.0 
TERMT(T.QE.TSTP) 

END  $ 

END  $ 

—TRANSLATION  TIME  -  O.BBT*** 
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Figure  9.  Computer  Program  Flowchart 
32 
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Figure  12.  Rotor  Speed  Versus  Resistance  (R^,) 
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'  Desired 
Horsepower 
(hp) 


Rotor  Speed  xIOOO  (rpm) 


Figure  13.  Rotor  Speed  Versus  HPD 
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open.  The  update  time  interval  of  0.1  second  is  within  the  control 
room's  test  computer  capabilities.  The  results  of  various  inputs  will 
be  discussed  in  the  next  chapter. 
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CHAPTER  IV 


COMPUTER  MODEL  RESULTS 

Model  Response 

The  computer  model  simulation  of  the  horsepower  Extraction  System 
was  run  with  step  and  profile  input  command.  The  step  input  of  100  hp 
command  shows  the  steady-state  accuracy.  The  profile  input  command 
establishes  the  speed  of  response  (reaction  time)  of  the  system  under 
transient  engine  operations.  The  results  of  the  model  simulation  are 
validated  against  test  data  for  a  rotor  speed  increase. 

Model  Validation 

Figure  14  compares  test  data  to  simulated  data  during  a  typical 
airstart  rotor  speed  increase  (360  rpm/sec).  The  agreement  between 
test  data  and  simulated  data  is  acceptable  within  the  range  of  avail¬ 
able  data  (20%  error  at  5000  rpm)  for  modeling  the  system.  The  simu¬ 
lated  data  plot  also  shows  the  maximum  horsepower  extraction  boundary 
for  each  rotor  speed  since  the  water  flow  is  pegged  at  its  limit  of  16 
gpm.  Only  the  rotor  speed  will  increase  the  amount  of  horsepower 
extracted. 

Step  Inputs 

The  response  to  a  100  hp,  constant  rotor  speed  set  to  10,200  rpm, 
step  input  command  with  a  one  percent  air  entrainment  is  shown  in 
Figure  15.  The  horsepower  settles  to  within  99%  of  the  command 
horsepower  level  after  2.0  seconds.  The  step  input  successfully 
demonstrates  that  accurate  steady-state  data  (99%  desired)  can  be 
obtained  within  a  few  seconds. 
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Rotor  Speed  x  1000  (RPM) 


Figure  14.  Comparison  Between  Simulated  and  Actual  Uaterbrake 
Data 
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Figure  15.  Response  to  a  100  HP  Step  Input  Command 
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Figure  16.  Response  to  a  100  HP  Step  Input  Command  with 
Reduced  Length 


The  model  response  to  a  step  Input  command  of  100  hp  with  a 
reduced  transmission  line  length  is  plotted  in  Figure  16.  The  reduced 
pipe  length  (1/10  of  full  length)  reaches  the  desired  horsepower 
approximately  0.25  seconds  faster  than  with  the  normal  line  length  but 
takes  about  three  seconds  longer  to  settle  99  percent  of  the  desired 
horsepower  (HPD).  This  short  line  length  does  not  dampen  out  the 
pressure  spike  from  the  valve.  The  waterbrake  capacitance  and 
inductance  which  was  not  Included  in  this  model  will  have  some  effect 
on  the  results  shown  in  Figure  16. 

Profile  Input 

The  response  to  profile  desired  horsepower  (HPD)  command  presented 
in  Figure  17  shows  the  dynamic  control  system  response.  The  engine 
rotor  speed  ramped  up  from  4000  rpm  to  10,000  rpm  at  a  rate  of  360 
rpm/sec.  The  profile  is  a  series  of  step  Inputs  and  a  ramp  input 
similar  to  aircraft  accessory  loads.  The  10  hp  step  command  increases 
at  4000  and  at  5000  rpm  are  both  reached  after  0.9  sec  overshooting  the 
desired  horsepower  by  approximately  3.5  hp.  The  step  input  of  20  hp  at 
6000  rpm  reaches  its  maximum  water  flow  horsepower  extraction  limit 
after  0.9  sec.  The  controller  easily  tracks  a  6  hp/sec  ramp  command 
from  7000  to  10,000  rpm  (at  360  rpm/sec).  The  control  system  satisfies 
the  secondary  goal  of  providing  relative  control  during  transient 
engine  operations  as  shown  in  Figure  17. 
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Rotor  Speed  xIOOO  (rpm) 


Figure  17.  Response  to  a  Profile  Input  Command 
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CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  waterbrake  horsepower  extraction  system  was  modeled  using  the 
lumped-parameter  technique  to  characterize  each  of  the  system  com¬ 
ponents.  A  program  was  developed  using  ASCL  to  simulate  the  system 
response  to  various  inputs.  The  results  show  that  both  of  the  original 
goals  were  satisfied;  to  obtain  accurate  steady-state  and  relatively 
accurate  and  quick  transient  horsepower  extraction  capabilities.  The 
computer  model  response  to  a  100  hp  step  input  shows  99^  of  the  desired 
horsepower  is  achieved  in  less  than  three  seconds  versus  manual  control 
that  takes  about  3-5  minutes  to  verify  steady-state  horsepower  level. 
The  profile  responses  show  the  system's  ability  to  track  a  series  of 
step  changes  and  a  ramp  input  ,  with  much  better  accuracy  than  an 
operator. 

It  is  recommended  that  Inductance  and  capacitance  be  Included  in 
the  waterbrake  model  to  improve  the  system's  dynamic  characteristics. 
The  inductance  and  capacitance  terms  can  be  quantified  by  measuring  the 
waterbrake  responses  to  known  inputs. 

Based  on  potential  cost  savings  in  turbine  engine  testing  time, 
the  system  could  recover  the  initial  cost  in  a  very  short  time.  In 
addition,  the  automated  horsepower  extraction  system  would  provide 
significant  improvement  in  controlling  horsepower  during  engine  air 
starts.  For  the  reason  stated  above. the  computer  controlled  extraction 
system  would  be  beneficial  to  implement  in  the  turbine  test  facilities. 
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